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The Mott scattering of high-energetic twisted electrons by atoms is investigated within the frame¬ 
work of the first Born approximation and Dirac’s relativistic equation. Special emphasis is placed 
on the angular distribution and longitudinal polarization of the scattered electrons. In order to 
evaluate these angular and polarization properties we consider two experimental setups in which 
the twisted electron beam collides with either a single well-localized atom or macroscopic atomic 
target. Detailed relativistic calculations have been performed for both setups and for the electrons 
with kinetic energy from 10 keV to 1000 keV. The results of these calculations indicate that the 
emission pattern and polarization of outgoing electrons differ significantly from the scattering of 
plane-wave electrons and can be very sensitive to the parameters of the incident twisted beam. In 
particular, it is shown that the angular- and polarization-sensitive Mott measurements may reveal 
valuable information about, both the transverse and longitudinal components of the linear momen¬ 
tum and the projection of the total angular momentum of twisted electron states. Thus, the Mott 
scattering emerges as a diagnostic tool for the relativistic vortex beams. 

PACS numbers: 03.65.Nk, 03.65.Pm, 34.80.Bm 


I. INTRODUCTION 

Beams of photons and charged particles that carry a 
non-zero projection of the orbital angular momentum 
(OAM) upon their propagation direction attract consid¬ 
erable attention in both, fundamental and applied re¬ 
search. During the last two decades, for example, a large 
number of experiments have been performed with optical 
twisted (or vortex) beams M- Moreover, very recently 
the production and manipulation of the twisted elec¬ 
trons have become feasible M- These vortex electrons, 
whose OAM projection may be as high as hm = 200 h Q, 
serve today as a valuable tool for probing the structural, 
electronic and magnetic properties with atomic resolu¬ 
tion. The engineering of the sub-nanometer sized vor¬ 
tices opens also new possibilities for studying the shell 
structure and dynamics of individual atoms &E3. 

Most applications of twisted electron beams to the 
study of (material) structure rely on the knowledge about 
the electron scattering by target atoms or ions. Accurate 
theoretical description of the basic scattering processes 
is required, therefore, which properly accounts for the 
details of the electron-atom interaction and the vortex 
structure of the incident beams. Moreover, such a the¬ 
oretical analysis has to be necessarily performed within 
the relativistic framework. This is due to the fact that 
the electron vortices are currently produced in transmis¬ 
sion electron microscopes 10 with the typical kinetic 
energy of 300 keV, corresponding to the velocity of about 
80 % of the speed of light. Even much higher energies up 
to 50 MeV are likely to be achieved in an experiment 
at the Thomas Jefferson National Accelerator Facility 


where the twisted electrons, created with the help of con¬ 
ventional holographic mask technique, are to be injected 
into the linear accelerator 0. The scattering of these 
MeV beams on solid targets, which will be observed to 
verify the vortex nature of the accelerated electrons, will 
proceed in an ultra-relativistic regime and definitely re¬ 
quires the proper theoretical treatment based on Dirac’s 
equation. 

Despite the definite need for the relativistic analysis 
of the twisted electron scattering by atomic targets, this 
field of research has barely started to develop. Up to 
the present, most of the collision studies with the vor¬ 
tex beams have been performed in the non-relativistic 
framework fl34l3 |. In contrast, the twisted electron so¬ 
lutions of the Dirac equation have been applied either in 
free space (TEI or in the presence of strong laser fields 
|16l - |l8j . but not in the scattering kinematics. Here we 
close this gap and investigate the Mott scattering of rel¬ 
ativistic twisted electrons. In our study, we pay special 
attention to the angular distribution and polarization of 
the outgoing electrons. In order to calculate these (angu¬ 
lar and polarization) properties we made use of the first 
Born approximation and the free-particle solutions of the 
Dirac equation. Moreover, we approximated the atomic 
potential by a sum of Yukawa terms; an approach that is 
justified for the description of the relativistic electron- 
atom collisions [19j. For such a potential, we evalu¬ 
ated in Section El the first-Born scattering amplitudes 
for both, initial plane-wave and twisted Bessel electrons. 
With the help of these amplitudes we were able to de¬ 
rive then the angular distribution and polarization of the 
outgoing electrons. For the incident vortex beam, more- 
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over, two experimental scenarios have been considered in 
which the target is (i) just a single well localized atom, 
or (ii) a macroscopic ensemble of randomly distributed 
atoms. The relativistic calculations have been performed 
for both scenarios as well as for the neutral hydrogen 
and iron targets, and incident electron energies in the 
range 10 < T e < 1000 keV. These calculations are dis¬ 
cussed in Section Hill and clearly indicate that the angular 
distribution and polarization of the scattered electrons 
are sensitive not only to the (ratio of) transverse and 
longitudinal linear momenta but also to the projection 
of the total angular momentum (TAM) of the incident 
vortex beam. Even though these results have been ob¬ 
tained within the first Born approximation, which has a 
rather limited accuracy for the description of (realistic) 
electron-atom collisions, they allowed us to elucidate the 
main effects arising in the relativistic Mott scattering of 
twisted beams. Our study, therefore, can be viewed as a 
starting point for more elaborate calculations including 
spin-orbit and spin-spin interactions as well as higher 
terms of the perturbative expansion. 

Relativistic units (H = 1 and c = 1) are used through¬ 
out the paper. 


II. THEORETICAL BACKGROUND 
A. Plane—wave electron scattering 

The (elastic) Mott scattering of plane-wave electrons 
has been worked out long ago and discussed in many 
textbooks. In the present paper, therefore, we will re¬ 
strict ourselves to a rather short compilation of the basic 
relations, which will be used later to treat the twisted 
electron scattering. 


1. Evaluation of the transition amplitude 

The total Mott cross section as well as the angular dis¬ 
tribution and polarization of outgoing (scattered) elec¬ 
trons are usually expressed in terms of the scattering 
amplitude. Within the framework of the first Born ap¬ 
proximation, this amplitude can be written as: 

h\'{p,p') = —^ J V’p' a' ( r )V( r ) V'pA (r) dr, (1) 

where V(r) is the electron-atom interaction potential. 
In Eq. ©> moreover, the wave-functions if p \(r) and 
ipp'y ( r ) describe a free electron in its initial and final 
states, characterized by the momenta p and p' and helic- 
ities A and X'. These functions are plane-wave solutions 
of the Dirac equation and read: 


where u p \ (and u p >y) is the Dirac bi-spinor [20l [2l| : 

( y/e ± m e u/ A )(n) \ 

_ I • (3) 

2\y/e — m e w^ (n) J 

In this expression, e = yj ml ± p 2 and n = p/p is the 
total energy and the propagation direction of an electron, 
and spinor w ^ (n) is the eigenfunction of the helicity 
operator A (n) = an/ 2, 

A(n)u/ A i (n) = w^ ( n ) = Xw^ ( n ), (4) 

with rr being the vector of Pauli matrices. As seen from 
this expression, the helicity A = ±1/2 is the projection of 
the electron spin onto its own direction of propagation. 
For the further evaluation of the scattering amplitude 

S one needs to find the explicit form of the spinor 
x \n). To achieve this goal, let us consider first the 
electron propagation along the quantization ( z -) axis of 
the overall system. In this case n = e z and the Eq. 6D 
simplifies to: 

^-w {a) (e z ) = s z w^(e z ) = au^ a) (e z ) , (5) 

thus indicating that w (xX \e z ) is just the standard Pauli 
spinor: 

«.(- 1,a '(e.) = U j (6) 



which describes the spin-up, a = ±1/2, and spin-down, 
a = —1/2, states along the z-axis. With the help of the 
solutions ([6]) we can construct now the spinor w^ (n) 
for the electron propagating in some arbitrary direc¬ 
tion n = (sin 9 cos <p, sin 0 sin ip, cos 9) with regard to the 
quantization axis: 

1/2 

w {x \n) = D ii(? , A 0 )ffl W (e z ) 

a——1/2 
1/2 

= E e- i ^d 1 J x 2 (0)w^(e z ). (7) 

er=-l/2 


Here = e latp is the Wigner D- 

function, and = cos (0/2) — 2ct<5 CTi _a sin (0/2), 

see Ref. 0. 

Making use of the expansion m and Eqs. © © we 
can finally derive the wave-function of the initial free 
electron state: 


V'pA(r) = e 


1/2 

ipr Q~ iCr,f d l Jy(0) 

cr=-l/2 

y/e ± m e (e z ) 


2Xy/e — m e (e z ) 


ip P x(r) = u p x e tpr , ifp'x (r) = u p >\> e lp r 


( 2 ) 


x 


( 8 ) 
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and similar expression for ?/ p ' A '(r). Beside these func¬ 
tions we have to define also the interaction potential V (r) 
in order to calculate the scattering amplitude (HJ). In 
the present study, we consider a spherically-symmetric 
Yukawa potential: 


where cos© = pp 1 /pp'. Apart from this scattering angle 
0, the cross section depends also on the initial and final- 
state helicities A and A'. If the incident electrons are 
unpolarized and the spin state of the scattered electrons 
remains unobserved, one derives from Eq. m-- 


Ze 2 

V(r) = -e -/ir , (9) 

r 

as approximation to the Coulomb field of the nucleus 
that is screened by atomic electrons, see Section ITTT1 for 
further details. By inserting this potential and the wave- 
functions ip p \{r) and if P 'y(r ), into Eq. ([TJ we find: 

/av(p,p) = q2 + 2 %, x , u pX 

2Ze 2 . . . 

= - 5 ——y {S0\\> +rn e d\.-x') 
q z + qr 

1/2 

x £ 0,(10) 
<t=-1/2 

where q = p — p’ is the momentum transfer, and in the 
last two lines we used the orthonormality of the Pauli 
spinors, (e z ) w^ a ^ (e z ) = 5 aa '■ In this expression, 
moreover, the angles (0, ip) and characterize the 

propagation directions of the incident and scattered elec¬ 
trons with regard to the overall coordinate system. 


2. Cross section and degree of polarization 


With the help of the amplitude (HOD we are ready now 
to investigate the properties of the plane-wave (relativis¬ 
tic) electrons scattered by the Yukawa potential. We 
start from the angle-differential cross section which is 
defined in the standard way as: 


d w 


2 p 

jz 


\h\'{p,p')\ 2 , 


( 11 ) 


where j z = ip P x{r)'y z ip p \(r) is the current density and j z 
is the Dirac matrix. For the initial-state wavefunction 
ip p \(r) as given by Eq. ©, the current density trivially 
simplifies to j z = 2p, and we find: 


d a 


(pi) 


^(B) = \fw(p,p')\ 
AZ 2 e 4 


(q 2 


P 2 ) 2 


(e 2 <^AA' + rn 2 S\ t -y) 


1/2 


r—-1/2 


2Z 2 e 4 

(q 2 +P 2 ) 2 


e 2 (l + cos©))(5aa' 


+ m 2 (l - cos0))(5 A _ A / , 


( 12 ) 


dcr o Pl) 

dfl' 


( 0 ) 


1 dcr AA / 
2^ dO 

AA' 

4^ 2 e 4 m 2 

{q 2 + p 2 ) 2 


( 0 ) 

(1 + A), 


(13) 


where the relativistic correction factor A is given by: 

A = --cos 2 (0/2) . (14) 

mi 


For p = 0 this formula reduces to the well known expres¬ 
sion for the Mott scattering by the Coulomb potential: 


dcTMott 

dD' 


(0') 


Z 2 e 4 

4 v 2 p 2 sin 4 (0'/2) 


(l — v 2 sin 2 (0'/2)) . 


(15) 

Here v = p/e and we assumed that the incident electrons 
propagate along the z-axis and, hence, 0 = 0 and 0 = 0'. 

The cross section m can be simplified also for the 
low-energy collisions, when T = e — m e <§( m e . In this 
case the correction factor A from Eq. (ED vanishes and 
the d4 pl) /df2' reads as: 


doo Pl) rn x = 4Z 2 e 4 m 2 
dfl' 1 ’ tf + p 2 ) 2 ’ 


(16) 


which corresponds to the non-relativistic scattering of 
spinless particles by an Yukawa potential. This simple 
formula can not be applied, however, to analyze present 
electron microscope (scattering) experiments in which 
the incident electrons have a kinetic energy of up to 
T = 300 keV. For such an energy, the correction factor 
ED increases to A « 1.5 cos 2 (0/2) which implies that 
both the total Mott cross section and the angular dis¬ 
tribution of scattered electrons can be strongly affected 
by the relativistic effects. This stresses the importance 
of the fully-relativistic analysis of the scattering of high 
energetic plane-wave (as well as twisted) electrons. 

The helicity-dependent cross section m can be ap¬ 
plied also to analyze the degree of longitudinal polar¬ 
ization of the scattered electrons P, which is also often 
referred to as the mean double helicity. For example, if 
the incident electron beam is completely longitudinally 
polarized, A = 1/2, we find: 


P (pl) (0) 


da. 


(pi) 

1/2,1/2 


— der 


d a 


(pl) 


dcr 


1/2,1/2 

R + cos 0 
1 + R cos 0 ’ 


(pi) 

1/2, -1/2 
(pi) 

1/2, -1/2 


(17) 
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where: 


R = 



(18) 


As seen from these expressions, in the ultra-relativistic 
regime R —> 1 and the helicity is conserved during the 
course of scattering, P(0) = 1. In contrast, for slow 
incident electrons R —> 0 and, hence, P = cos© which 
corresponds to the projection of the initial electron spin 
onto the final momentum p'. 


B. Twisted electron scattering 

After having briefly recalled the basic relations used to 
describe the Mott scattering of plane-wave electrons by 
the central Yukawa potential, we are ready to consider 
the twisted electron beam. Similar to before, we shall 
start from the derivation of the scattering amplitude and 
discuss its main properties. 


= nr 


dpp 


_^ e ^ p -ip x b /AA/(p;p , )(21 ) 

Z7T 


where fx\'{PiP r ) is the “standard” plane-wave matrix 
element (fTUl) . In Eq. (1211) . moreover, we have introduced 
the exponential factor exp(— ip±b) to specify the lateral 
position of the scatterer atom with regard to the cen¬ 
tral ( 2 -) axis of the incident electron beam, and where 
b = (b x ,b y , 0) is the impact parameter. The introduc¬ 
tion of this factor reflects the fact that, in contrast to the 
plane-wave case, the Bessel beam has a complex spa¬ 
tial structure in the plane perpendicular to the z-axis 
01 - 

In order to further evaluate the F^ (p,p', b) we have 
to insert the plane-wave amplitude m into Eq. ED- 
After performing some simple algebraic manipulations, 
we finally obtain: 


Fix\p,p', b ) = 


2Ze 



+ m e Sx,-y) 


1. “Twisted” transition amplitude 


In order to derive the amplitude for the Mott scatter¬ 
ing of twisted electrons we have to return to the general 
expression ©• In this formula, the initial-state wave- 
function ijj p x(r) should be modified to represent a twisted 
state. Here we assume that the incident twisted elec¬ 
trons propagate along the quantization (z-) axis and have 
well-defined values of (i) the longitudinal linear momen¬ 
tum p z , (ii) the modulus of the transverse momentum 
|p_i_| = x, and (iii) the half-integer projection of the to¬ 
tal angular momentum, J z = m. Such a Bessel state 
has, moreover, the definite energy e = -^/x 2 + p 2 + ro 2 
and helicity A, and is described by the wave-function: 


1p>cmp z A ) 


d 2 RL 

(2tt) 2 


axm(p±)u p \ e lpr , 


(19) 


see Refs. [13,113,113 and Appendix [A] for further details. 

As seen from this expression, the ip^mp z x(r) can be con¬ 
sidered as a coherent superposition of the plane-waves 
Up\d lpr , weighted with the amplitude: 


a„m(p±) = H) m e imVp £ 5 (|p ± | - x) . (20) 

V x 

The linear momenta of these plane-wave components, 
p = (p±,p z ) = {xcostfip, xsin(pp,p z ), form the surface 
of a cone with the opening angle 9 p = arctan (x/p z ). 

Employing the initial-state wave-function (1191) we find 
the amplitude for the Mott scattering of the Bessel elec¬ 
trons: 

F^{p,p', b ) = -±J^Ar)V(r)i>« mp Ar)dr 


d 2 p± 

H) 2 


(p ±) e WA - h\>(p,p') 


1/2 

xe-V £ d 1 J x 2 (e p )d 1 J 2 (d')I m . <7 (a,l3,b),(22) 

cr=-l/2 

where the function / m _ CT (o:, /3, b) reads as: 


2tt 


I n (a,f3,b) = J 


d 4> e 


i n4>—ixb cos(</>+<^ —(fib) 


27T 


a — p cos < 


(23) 


In this expression, ipi, is the azimuthal angle of the 
impact-parameter vector b, and we have introduced the 
short-hand notations: 


a = 2p 2 - 2p.p' z + p 2 

= 2p 2 {l — cos d p cos 9') + p 2 , (24) 

/3 = 2xp' ± = 2p 2 sin 0 p sin 0', (25) 

where both a and [3 are positive. 

Eqs. (E^l) ED represent the amplitude for the Mott 
scattering of twisted (Bessel) electrons by the Yukawa 
potential ED, whose center is shifted by b with respect 
the z-axis. In the following Section we will apply this 
amplitude to study the angular distribution of the scat¬ 
tered electrons. 


2. Impact-parameter dependent angular distribution 

In contrast to the plane-wave case CD, the evaluation 
of the cross section for the scattering of a vortex beam 
by a single well localized atom is not a simple task. It 
requires the re-definition of the incident flux j z . which 
is now a function of the impact parameter b and, more 
generally, the concept of the cross section. The deriva¬ 
tion of such an impact- parameter-dependent Mott cross 
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section is out of scope of the present study and will be 
discussed separately [25!]. Here, instead, we will focus on 
the angular distribution of the scattered electrons: 


w® (*V; b ’ m ) = * F^( PlP r ,b) 


= a r- 


■ x(e 2 5\\> +m 2 e 6\-\') 


1/2 


r=-l/2 


(26) 


where the pre-factor A f is defined by the normalization 
condition J (6', <p' b,m) dfT = 1. 

From the discussion above and Eq. (l26l) . it is apparent 
that the angular distribution (9', ip'; b,m ) of the 

outgoing electrons depends on the position of the Yukawa 
scatterer within the incident wave front. Even though the 
detailed study of such a 6-dependence requires numerical 
computations and will be presented in Section m here 
we consider the special case of a central collision, 6 = 0. 
For this scenario, one can evaluate the scattering ampli¬ 
tude (l22l) and, hence, the W’jf'y (O',<p'] b,m) analytically. 
In particular, we re-write the integral (1231) . which enters 
in F w(P,P', b ), as: 


J„(a,/3,6 = 0) 


n<f) 


d <j> 

2n a — /3 cos <f> 


A n \ 


iir/3 / (zi - z)(z - z 2 ) 


dz, (27) 


where in the second line we introduced the complex vari¬ 
able z = e 1 ^. The integration contour in Eq. (1271) is the 
unit circle and the positions of the poles are given by 
Zip = (a ± \Ja 2 — /3 2 ^ / fd with z\ > 1 and z 2 < 1. By 

noting that the integrand in Eq. m has a single simple 
pole inside the unit circle, and taking its residue at z 2 . 
we find: 


I n (a,P,b = 0) 




M 


\Ja 2 — 13 2 \ a + \/a 2 — / 3 2 


(28) 

With the help of this expression one can easily evaluate 
the angular distribution (l26l) of the scattered electrons 
for zero impact parameter 6 = 0. Of special interest 
here is the behaviour of (8' b = 0,m) at small 

scattering angles with respect to the incident beam axis. 
If we perform the Taylor expansion of Eq. (l26l) for 8' —> 0 
we obtain: 


WS?(0', if/-, 6 = 0, to) cx (0') 2|m V| , (29) 


which indicates that the angular distribution of the scat¬ 
tered electrons vanishes for the forward emission if the 


incident beam is twisted, 6 = 0 and to ^ A'. This be¬ 
haviour has been predicted recently also within the non- 
relativistic framework ITol , and is strongly different from 
what is expected for the “plane-wave” angle-differential 
cross section m which is maximal for 8' = 0. 

The dip (l29l) in the electron angular distribution for 
8' = 0 disappears with the increase of the impact param¬ 
eter. Indeed, by making use of Eq. (177|) we find that the 
angular distribution: 


<v(<?' = (W;6,to) = U L 

’ 7r(4/rsm (8 p /2) + p z ) z 

x [e 2 (l + cos 9 p ) 6 X \' + TOg(l - cos 9 p ) 5 x ,_a'] , (30) 


1 Z 2 e 4 x [Jm-\'{>cb)Y 


does not vanish for the forward scattering and 6^0. For 
small impact parameters, b <C l/*c, this expression pre¬ 
dicts that the electron emission quickly increases with the 
impact parameter, (9' = 0,<p'; b,m) oc b 2 


3. Averaging over the impact parameter 


Until now, we have discussed the evaluation of the an¬ 
gular distribution of scattered electrons for the case when 
the incident twisted beam collides with a single and well 
localized Yukawa potential. This scenario, however, can 
be hardly realized in nowadays Mott scattering experi¬ 
ments in which thin foils are usually used as a target. We 
can describe the (solid-state) target by an ensemble of 
Yukawa scatterers that are randomly and uniformly dis¬ 
tributed over the transverse extent of the incident beam. 
If one neglects collective and multiple scattering effects, 
the differential cross section for the Mott scattering of 
twisted electrons by such a “foil” can be introduced in 
a way similar to that of the standard plane-wave case 
m■ That is, by averaging the (square of the) transition 
amplitude (1211) over the impact parameter 6 we obtain: 



where J z is the total incident current: 


Jz = J Jz(b)d 2 b 

= / ^>cmp z \{r - b)^ z i!) >tmpz x{r - 6)d 2 6, (32) 


and ipxmp z \{ r ~ 6) is the wave-function of the electron 
vortex beam whose axis is shifted by 6 with respect to 
the overall z-axis. 

In order to further evaluate the Mott scattering cross 
section we should perform the integration over the impact 
parameter 6 in Eqs. (1311) and (1371) . Assuming that the 
radius R describes the characteristic transverse size of 
the Bessel beam and R /?> l/x, the total current J z can 
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be approximated as: 


\Jz — I 'l/j?cmp z \( t>) 'lp*cmp z \{ 0)(\ I) 


— < ^‘^ Z j 


d RL O R / Q Q \ 

= 2 p z —, (33) 


(2tt) 


where in the second line we have employed Eq. (76a) from 
[26| to evaluate the square of the 8 function, 


\S(\p ± \-x)\ 2 =-S(\pa_\- x). (34) 

7 r 

By making use of this expression and of Eq. m one can 
also compute the integral of the square of the transition 
amplitude: 


/ 



R 


7r 


2tt 


/ 


\f\\'(p,p' 


2 d ifip 

2t r ’ 


0 


(35) 


and, hence, the angle-differential cross section: 


j ( tw ) 
dcr AA' 

dll' 


(o', e p ) 


i 

COS dp 


J \h\'(p,p ')\ 2 

0 


(36) 


To carry out the remaining integration over ip p we in¬ 
sert here the explicit form of the plane-wave transition 
amplitude (fTOl) : 


du 


(tw) 


AA 


dO 


d v) = 


4Z 2 e 4 

COS 0 V 


(e 2 5\x' + m 2 e 8\-\') 


\d 1 ^{0p)d 1 J^0 l )d 1 J,l{0 p )d 1 J,l\e') 


J e i ( CT - CT ')(^-^') dipp 

J ( q 2 + /r 2 ) 2 2 tt 

0 


(37) 


and observe that the integral in the last line of this ex¬ 
pression is given by: 


2tt 


0 


e i(<r-<r')(<p p -<p') dip p 

(q 2 + /i 2 ) 2 27 t 


ot 8(7(7' T fd8(7^ — (7' (OO) 

(a 2 - /3 2 ) 3 / 2 ’ 1 ' 


where the function I a i- a (a,l 3,0) and coefficients a and 
are defined by Eqs. (1251) . (l24l) and (1251) . respectively. 
By inserting the integral (1551) into Eq. (EH we derive 
the final expression for the angle-differential cross sec¬ 
tion of the Mott scattering of the twisted electrons by a 


macroscopic target: 


d a 


(tw) 


AA 


dft 


H*', Op) = 


2„4 


4Z 2 e 

COS 0 n 


(e 2 8xy + m 2 8x-x') 


x E 4a(*p) &')rtl(° P ) 


cr 877' -f- {8 8 ( 7 ^—c 

(a 2 - /3 2 ) 3 / 2 


(39) 


where the polar angle O' of the outgoing electrons is de¬ 
fined with regard to the z-axis. As seen from Eq. m, 
the duS^/dO' is insensitive to the TAM projection m 
of the (initial) Bessel beam but depends on its opening 
angle 0 p . If 0 p = 0 and, hence, a = 4p 2 sin 2 (0'/2) + p 2 , 

(0 = 0, and d aX (0 P ) = S a x, this cross section expectedly 
coincides with the plane-wave result (fl2l) where 0 = 0'. 

The angle-differential cross section (1551) still depends 
on the helicities of the initial and final-state electrons, A 
and A'. It can be used, therefore, to evaluate the angular 
and polarization properties of the scattered electrons for 
every possible experimental setup. In Section El for 
example, we will discuss the angular distribution of the 
scattered electrons under the assumption that their spin 
state remains unobserved. For this case, the Mott cross 
section: 


Hrr (tw) rU tw 

V«Mp> = e^'.'w 


dft' 


(40) 


2 Z 2 e 4 


( a 2 — /3 2 ) 3 / 2 cos 0. 


V L 


(e 2 + m 2 ) 


+p 2 (a cos(0 p /2) cos(0'/2) + (3 sin(0 p /2) sin(0'/2)) 


is insensitive to the helicity of the incoming beam and, 
within the limit p 2 -C m 2 , reduces to the non-relativistic 
expression: 

dfi' e,er ' cos dp {a 2 -P) 3 ' 2 ' ( ’ 

which describes the scattering of spinless particles by the 
Yukawa potential [25]. 

Similar to the plane-wave case, we can employ Eq. ESD 
also to evaluate the degree of longitudinal polarization of 
the outgoing electrons. This degree p( tw )(0',0 p ) is given 
by the first line of Eq. m where the plane-wave cross 
sections are substituted by dcrS^/dO'. For the sake of 
shortness, we will not present here the final expression 
for the p( tw )(0',0 p ) but discuss in Section m the nu¬ 
merical predictions for the polarization of the final-state 
electrons. 


III. SCATTERING BY ATOMIC POTENTIALS 

Up to the present we have discussed the Mott scatter¬ 
ing of plane-wave and twisted electrons by the Yukawa 
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FIG. 1: (Color online) The angular distribution (1471) of the electrons scattered by a single hydrogen atom. Calculations have 
been performed for the incident Bessel electron beam with the kinetic energy 10 keV, helicity A = 1/2, opening angle 6 P = 20 deg 
(top panels) and 9 P = 40 deg (bottom panels), and the projection of the TAM m = 1/2 (blue dashed line), m = 3/2 (green 
dash-dotted line) and m = 5/2 (red dotted line). Moreover, the target atom is assumed to be displaced by distances 6 = 0 nm, 
b = 0.01 nm, and 6=1 nm with respect to the central beam axis and under the azimuthal angle tpi = = 0 deg. Results of 

the calculations are compared with the prediction obtained for the incident plane-wave electrons (solid black curve). 


potential V(r) given by Eq. (0. In atomic collision 
studies, the V(r) is often used to approximate (realis¬ 
tic) electron-atom interactions. For example, the elec¬ 
trostatic potential of a neutral atom can be written as a 
sum of three Yukawa terms: 

7p 2 3 

U at (r) = — (42) 

i=1 

where the parameters A and //,; are determined by a 
fitting to the results of the Dirac-Hartree-Fock-Slater 
(DHFS) self-consistent calculations jl9|, j23|. It is rather 
straightforward to generalize all the results obtained in 
the previous Section to the case of such an “atomic” po¬ 
tential. Indeed, the scattering amplitudes can be written 
in this case as: 


fxx'(p,p') = 2 Ze 2 (£ 


A, 


\i=1 


q 2 + p 2 


{ebxx' + m e 6\-\f) 


1/2 


£ e^'-^y>)dy>'), (43) 

<7= —1/2 


for the plane-wave, and: 

Fx>/{p,p', b ) = 2Ze 2 i~ m (e5\x+m e 5 X -x') 


1/2 

xe <v £ d^iejdl^e') 

a—— 1/2 
3 

x £ A i I m _ <7 (a i ,P, b) , (44) 

i= 1 


for the twisted-wave electrons, where = 2p 2 (l — 
coscos#') + /if. With the help of these amplitudes 
one can again evaluate analytically the angular distri¬ 
bution and the (degree of) longitudinal polarization of 
scattered electrons both for the plane- and twisted-wave 
cases. For example, by inserting the expression (H51) into 
Eq. (1M1) we find the angle-differential cross section for 
the scattering of Bessel electron beam on a macroscopic 
atomic target: 


df>' 


(#', e p ) 


4Z 2 e 4 
cos 6p 


X 


{e 2 5xw + mlSx-y) 











x E4 a (*p) 4v(*')4£(«w4'v(0O 

Target 

Ai 

Hi 

^4.2 

H2 

M3 

era' 

3 

H (Z=l) 

-184.39 

2.0027 

185.39 

1.997 

0 

x ^ ^ Ai Ak X,j'(ai , a/-, /3), (45) 

i.k=1 

Fe (Z=26) 

0.0512 

31.825 

0.6995 

3.7716 

1.1606 


where the integral T a >- G (pLi, aj3) reads as: 

2tt 


? /^) — ^ 


pincf) 


d(j) 


0 

2tt 


(a* — f3 cos 4>)(ak — /3cos0) 27r 


Q m0 


pin0 


cx.k — oti J \ai — f3 cos0 a/c — j3 cos </>/ 27r 
o 

_ J ra (ctj,/3,0) - I n (a k ,l3, 0) 
and I n (a,j3, 0) is given by Eq. 


d</> 

?7T 

(46) 


IV. RESULTS AND DISCUSSION 


In the previous Sections we have derived the Mott cross 
section and angular distribution of the outgoing elec¬ 
trons for two cases in which the incident Bessel beam 
collides with either (i) a single atom, or (ii) a macro¬ 
scopic target consisting of randomly distributed atoms. 
Even though the second scenario can be more easily re¬ 
alized experimentally and, hence, is of definite practical 
interest, let us first discuss the interaction of the twisted 
wave with a well-localized (single) atom. In Fig. [IJ for 
example, we display the angular distribution of electrons 
scattered by the neutral hydrogen atom and observed by 
a polarization insensitive detector: 

Rf w) 09 b,m) = J2 ( 0, > <p' = 0; b, m) . (47) 

A' 


Here, W^{8'^'- b,m) = A f 


F \w {P,P', b ) and the 


amplitude F is given by Eq. (|44| with the parameters 
Hi and Ai listed in Table [T| Calculations have been per¬ 
formed for the incident Bessel beam © with the energy 
10 keV, helicity A = 1/2, opening angles 8 p = 20 deg 
(top panels) and 9 p = 40 deg (bottom panels), and three 
projections of the total angular momentum: m = 1/2 
(dashed line), m = 3/2 (dash-dotted line) and m = 5/2 
(dotted line). We assume that the hydrogen atom is 
placed at the azimuthal angle (fib — <p' = 0 deg and dis¬ 
tances 6 = 0, 0.01 and 1.0 nanometers with respect to 
the beam axis. As seen from the figure, the angular dis¬ 
tribution (1471) is very sensitive to the parameters of the 
incident twisted wave as well as to the position of the tar¬ 
get atom. For 6 = 0, for example, the forward electron 
scattering is possible only if the TAM projection of the 
incident beam is m = 1/2. For m = 3/2 and 5/2, in con¬ 
trast, the W^ w \9'\ b,m) vanishes identically if 9' —> 0 
as expected from Eq. (1291) . If the target atom is shifted 


TABLE I: Parameters of the effective potential ([42jl used in 
the present calculations. The parameter A 3 = 1 — A\ — A 2 . 
Results are from Ref. l23ll . 


from the center of the incident beam, the angular distri¬ 
bution (1471) exhibits qualitatively similar 9' behaviour for 
all values of m. While a forward scattering becomes al¬ 
lowed for 6 > 0, c.f. Eq. (13U1) . most electrons are emitted 
under the polar angle 9' ~ 8 p where the 9 P is the open¬ 
ing angle of the incident twisted beam. One understands 
such a 8' behaviour of the angular distribution by noting 
that the (initial) Bessel state can be seen as a coherent 
superposition of plane waves lying on a momentum cone 
surface with the opening angle 9 pi see Eq. (fl9l) . The 
scattering of a plane-wave component e lpr by the atom 
results in the emission of the outgoing electron along the 
linear momentum p 1 which is tilted by the angle 9 p with 
respect to the overall z-axis. 

Up to the present, we have discussed the calculations 
for the scattering of the twisted electrons by the sin¬ 
gle well-localized hydrogen atom. In more experimen¬ 
tally realistic scenario, the incident Bessel beam collides 
with a macroscopic target. By making use of Eq. (l45l) 
and parameters given in Table [T] one can investigate the 
angle-differential Mott cross section for this (realistic) 
case. For example, Fig. [2] displays the results obtained for 
the iron target and kinetic electron energies in the range 
10 < T e < 1000 keV. Similar to before, we assumed that 
the polarization state of the scattered electrons remains 
unobserved and, hence, summed over the helicity A': 


d 4 tw) 

dfy 


(O', 0p) 


E 

A' 


j (tw) 
dg ~AA' 

MV 


(0',e p ). 


(48) 


As mentioned already in Section [TTJ the averaging over 
the impact parameters of the target atoms leads to the 
fact that this cross section do- dwd /df V is insensitive to the 
projection m of the TAM and depends only on the open¬ 
ing angle of the incident Bessel beam 9 p . To explore the 
^-dependence, here we present results for three opening 
angles: 9 p = 5 deg (dashed line), 20 deg (dash-dotted 
line) and 40 deg (dotted line). As seen from the figure, 
the d<7 d '' d /dU / is again peaked at the angle 9' = 9 P and 
the effect becomes more pronounced with the increase of 
the incoming electron energy. This is in contrast to cal¬ 
culations for the incident plane-wave beam which pre¬ 
dict the predominant forward emission, c.f. solid line 
in Fig. [2] Therefore, the most remarkable distinction be¬ 
tween the angle-differential cross sections for the scatter¬ 
ing of twisted and plane-wave electrons can be observed 
for large values of the angle 9 p . 
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FIG. 2: (Color online) The angle-differential cross section 
(HO) for the Mott scattering of electrons by the macroscopic 
iron target. Computations have been performed for the 
incident Bessel beam with the helicity A = 1/2, energies 
T e = 10 keV (upper panel), 100 keV (middle panel) and 
1000 keV (bottom panel), and the opening angles 0 V = 5 deg 
(blue dashed line), 20 deg (green dash-dotted line) and 40 deg 
(maroon dash-double-dotted line). Results of the calcula¬ 
tions are compared with the prediction obtained for the inci¬ 
dent plane-wave electrons (solid black curve). 

We can use Eq. (l45l) to explore not only the angle- 
differential Mott cross section but also the degree of lon¬ 
gitudinal polarization of scattered electrons. If we insert 
do'AA' 5 /d^2 / (6 >/ ; into Eq. (flTl) we find the polarization 
of scattered electrons P^ tvJ \9]9 p ) for the interaction of 
the Bessel beam with a macroscopic atomic target. We 
have computed this degree of polarization again for the 
iron target and for the same beam parameters as in Fig. [2] 
Results of our calculations are presented in Fig. [3] and 
indicate that the p 6 w )( 0 - dp) can be rather sensitive to 
the opening angle 9 p of the Bessel beam. In particu¬ 
lar, while for the scattering of the plane-wave electrons 
with helicity A = 1/2 the polarization decrease mono- 
tonically from p( p b = l for 9' = 0 deg to p( p b = — 1 
for 9' = 180 deg, the p 6 w )( 0 ; 0 p ) is maximal at 9' = 9 p 
and falls down for the forward and backward emission. 
The influence of the opening angle 9 p on the longitudinal 
polarization of outgoing electrons can be observed most 



FIG. 3: (Color online) The degree of longitudinal polarization 
of electrons scattered by the macroscopic iron target. Calcu¬ 
lations have been performed for the same set of parameters 
as in Fig. [2] 


easily for low collision energies and relatively small scat¬ 
tering angles, 9' < 60 deg. In this parameter range both 
the Mott cross section and the variation of p( tw ) ( 9; 9 p ) 
with 9 p are large. For example, for T e = 10 keV and 
the observation angle 9' = 40 deg, the degree of polar¬ 
ization increases from 76 % to almost 100 % if instead 
of the plane-wave- the twisted-wave electrons with the 
opening angle 9 p = 40 deg collide with the target. 

In the calculations above we have always assumed that 
the incident electrons are prepared in a pure quantum- 
mechanical state with the well-defined projection of the 
total angular momentum m on the propagation (z-) axis. 
During the recent years, however, a number of experi¬ 
ments have been performed and/or proposed with the 
coherent superposition of two (or even more) states with 
different TAM projections [2 [H[. For two Bessel 

beams with the same helicity A and the same beam axis, 
such a superposition is described by the wave-function: 

V’^tw)(r) = ciVw lP *A(r) + c 2 Vw 2P , x(r) , (49) 

where ipxm n p z \( r ) is given by Eq. (fl9l) and c n coefficients 






















read as: 


c n = \cn\e loln , |ci| 2 + |c 2 | 2 = 1. (50) 


With the help of this expression and Eq. (1M1) we can find 
the angle-differential cross section 


dO' 


(O', 0 P ) 


(51) 


i 

cos 6p 


2tt 

J |/aa'(p,p ')| 2 G(ip P ,Am,Aa)^ 

o 


for the Mott scattering of the (superposition of) Bessel 
beams by a macroscopic atomic target. In contrast to 
Eq. (1361) . derived for a single incident beam, this cross 
section depends on the function 


G(<p p , Am, Act) = 1+2 |ci||c2| cos (A m(<p P - 7t/ 2) + Act)) , 

(52) 

and, hence, on the differences of the TAM projections, 
Am = to 2 — mi, and phases, Act = ct 2 — a\, of two 
twisted states. Such a Am- (as well as Act- ) depen¬ 
dence translates directly into the angular and polariza¬ 
tion properties of scattered electrons. In Fig. [2 for ex¬ 
ample, we display the cross section da^ tw '/df l'(O', 0 p ) = 

/ diY (O', 6 p ) and the degree of longitudinal po- 
A' 

larization p( 2tw ) (0; Q p ) for the scattering of 300 keV elec¬ 
trons by the iron target. Calculations have been per¬ 
formed for the opening angle 0 p = 40 deg, initial helicity 
A = 1/2, and for Am = 1 (solid line), Am = 2 (dashed 
line) and Am = 3 (dotted line). As seen from the figure, 
both the cross sections and the polarization are strongly 
affected by the variation of the Am. For example, the 
dcr^ 2tw ) /dSl'(0', Op) decreases by about factor of 4 if ob¬ 
served at the angle O' = 0 p = 40 deg and the difference of 
the TAM projections changes from Am = 1 to Am = 3. 
Such a remarkable effect can be easily observed by mod¬ 
ern electron detectors and may provide useful informa¬ 
tion about the scattering of (superpositions of) Bessel 
electron beams. 


V. SUMMARY AND OUTLOOK 

The first Born approximation and Dirac’s relativistic 
theory have been applied to explore the Mott scatter¬ 
ing of high-energetic twisted electrons by atoms. In our 
study, we focused especially on the angular distribution 
and longitudinal polarization of outgoing electrons. To 
derive these—angular and polarization—properties we 
have approximated the atomic potential by (a sum of) 
Yukawa terms and considered two different “experimen¬ 
tal” setups. In these setups, the incident Bessel beam 
collides with either (i) a single well localized atom, or (ii) 
a macroscopic target, consisting of randomly distributed 
atoms. In the first case, we found that the angular distri¬ 
bution of the outgoing electrons depends not only on the 




Polar scattering angle 0’ (deg) 


FIG. 4: (Color online) The angle-differential cross section 
(upper panel) and the degree of longitudinal polarization 
(bottom panel) of the electrons scattered by an iron target. 
Calculations have been performed for a coherent superposi¬ 
tion of two equally weighted Bessel beams with the difference 
of the TAM projections Am = 1 (black solid line), Am = 2 
(blue dashed line) and Am = 3 (red dotted line). The beams 
are prepared, moreover, in the state with the relative phase 
Aa = 60 deg, opening angle 9 P = 40 deg, helicity A = 1/2, 
and kinetic energy T e = 300 keV. 


kinematic parameters of the twisted wave, such as the ra¬ 
tio of the transverse x to longitudinal p z momenta, but 
also on the (projection of the) total angular momentum 
to. The most pronounced effect of the TAM can be ob¬ 
served if the (single) target atom is placed close to the 
centre of the twisted wave front. 

In contrast to the collision with a well localized atom, 
the angle-differential cross section for the Mott scattering 
of a single Bessel wave by a macroscopic target appears 
to be independent of the topological charge to. Never¬ 
theless, the (angular and polarization) properties of the 
outgoing electrons are still sensitive for this case to the 
beam’s opening angle, 0 p = arctan(>«:/p 2 ). In particular, 
the angle 0 p determines the direction of the predominant 
emission of the outgoing electrons. Such an emission pat¬ 
tern, that is peaked at O' = 0 p , provides a clear signature 
of the Mott scattering of twisted electrons and can be 
easily observed experimentally. 

For a macroscopic target we have also analyzed the 
scenario in which the incident electron beam is prepared 
as a coherent superposition of two Bessel states. The 
use of such a superimposed beam helps to restore the 
sensitivity of the Mott process on the projections m of 
the TAM. In particular, our calculations clearly indicate 
that the angular and polarization properties of the scat- 
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tered electrons strongly depend on the difference mi— m 2 . 
This Am-dependence can be observed most clearly at 
the scattering angle 9' = 9 p at which both the angle- 
differential Mott cross section and the longitudinal po¬ 
larization of scattered electrons are maximal. 

The present theoretical study has been performed in 
the (relativistic) first Born approximation. While the 
accuracy of this approach is limited, especially for rela¬ 
tively low collision energies and backward scattering an¬ 
gles, it helps to elucidate the main features of the elastic 
scattering of twisted electrons. Moreover, based on the 
developed theory, one can perform a more elaborate anal¬ 
ysis of the Mott process, including the spin interaction 
effects and higher perturbation terms. Such an analysis 
is currently underway and its results will be published 
elsewhere. 


Acknowledgements 

We are grateful to G. Kotkin and D. Karlovets for 
useful discussions. V. G. S. acknowledges support from 
RFBR (via grant No. 13-02-00695) and MES (Russia). 


Appendix A: Twisted electron wave—function 

Apart from this study, the wave-function of relativis¬ 
tic twisted electrons has been utilized in several recent 
works [HI, Ojl, UU . Since the explicit form of this function 
slightly differs in each paper, let us discuss and compare 
here the various representations of corresponding wave- 
functions. We start from our definition (fl9l) which, upon 
substitution of the plane-wave solution © and the am¬ 
plitude (f 20 ]l . reads as: 

ipxmpAr) =Y,j ^7\/^e i(, ’" +xr - LCOB( ' Pp "' Pr)) 


x (- i) m eHm-^Vp d y A 2 (d p ) U {r7) (e, A), (Al) 

where we performed trivial integration over p± and, for 
the sake of shortness, introduced the notation: 

/ \Z £ + m eW W (e*) \ 

f/ W (M)= _ . (A2) 

\ 2Xy/e - m e w^\e z ) ) 

This bi-spinor (e, A) is obviously the eigensolution 
of the spin operator 



with the eigenvalues a = ± 1 / 2 . 

By using the well-known integral relation 


we can perform the integration over the azimuthal angle 
ip p in Eq. (IA1I) and obtain the twisted wave-function: 

VWA(r) = J^e i ^Y,r aei{rn ~ aWrd lx^ 

” a 

x J m _ CT (^rj_) £/ (< 7 ) (e, A), (A5) 

in cylindrical coordinates (rj_, <p r % z )- The normalization 
condition for these states is 


J il>xmp z \{*)4>x'm'p^\'(r)d i r = S(x - V) 5 mm ' (A6) 

x 2ttS(p z - p' z ) ■ 2m e 5w . 

Each summand in the right- hand side of Eq. (IA5I) is the 
eigen-function of the operators L~ = —id/dip r , £ z and 
J z = Lz ± with eigen-values to — <7, a and to, respec¬ 
tively. However, due to the summation over a = ±1/2 in 
Eq. (IA5I) . the 4’xm Pz \(r) does not possess definite (z- ) 
projections of the spin and orbital angular momentum 
separately. It is the eigen-state of the total angular 
momentum operator, J z 4’xm Pz \(r) = with 

half-integer m. Therefore, the Bessel vortex state (ESI) 
(or (EH)) of free relativistic electrons includes intrinsi¬ 
cally spin-orbit interaction fl5| . 

The coordinate representation of the wave-function 
dm can be easily generalized to the case when the vor¬ 
tex electron beam is shifted by the vector b = ( b x , b yi 0 ) 
with respect to the quantization z-axis. Inserting the 
translation factor exp (—ipb) into Eq. (fl9l) and perform¬ 
ing integration over the p ±, we find: 

4>xm Pz \(r,b) = JZ^Y, i ~ CTei{rn ~ CT) ' Pp d 1 Jx^r) 

’ U 

x J m _ a (x P± )U^\e, A), (A7) 

where p = r± — b = (p±,<p p , 0 ) is the displacement vector 
in the x — y plane written in cylindrical coordinates. An 
alternative representation of the (displaced) Bessel elec¬ 
tron wave-function was prop osed in the non-relativistic 
framework in Refs. [U [29j . It relies on the addition 
theorem for Bessel functions and, although being math¬ 
ematically equivalent to Eq. JA3, is less convenient for 
practical purposes and does not allow a transparent phys¬ 
ical interpretation. 

By making use of Eqs. (IA5I) and dm , it is easy to 
investigate the behaviour of the Bessel electron wave- 
function for vanishing values of the transverse momen¬ 
tum, x —> 0, and, hence, the opening angle 9 p —> 0. In 

this limit d,J x (9 p ) = 5 a \ and J m _ cr (xp_ l) = Sm-< 7 , 0 , and, 
hence, 




xmp z X 


( r , b) 


X 


>«r—>-0 



Ve ipr ^=A, (AS) 


1 /' 27r 

/ e imp+iz cos v ^ 

27T J q 


i n Jn(z ), 


(A4) 


where the u p \ is the standard Dirac bi-spinor (|5|). and 
p = (0,0, p z ). As can be seen from this expression, the 
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Bessel electron wave function recovers—for x. —> 0—the 
standard solution for a plane wave that propagates along 
the z-axis with the helicity A = m = ±1/2. 

As mentioned already above, some another representa¬ 
tion of the wave-function ifxrnp z x(r) has been proposed 
recently by Karlovets [Tb] which reads as: 

^m' Pz x(r) = 

V a 

x Jm'-A.xr^) U (a \e, A), (A9) 

which differs from Eq. (IA5I) by (i) the overall pre-factor 
exp(iAyv) and (ii) the missing term i~ a under the sum¬ 
mation sign. The first factor leads to the fact that the 
ipxm'pzxir) is the eigen-function of the total angular mo¬ 
mentum operator J z with the eigen-value m! ± A, where 
m! is the integer. This factor comes just from the dif¬ 
ferent convention in the definition of the free-electron 
spinor, w^ x \n) = exp(i\tp)w( x \n), where u/ A )(n) is 
given by 0 , and does not affect—upon relabeling of the 
quantum numbers—the Mott scattering cross section. In 
contrast, the term i~ a which defines the relative phase of 


two components of the wave-function (IA9D with spin pro¬ 
jections a = ±1/2, was mistakenly missed in Ref. [161 ]. 
Without including this factor, the ip 3{ m'p z x(.r) is not a 
solution of the Dirac equation with a definite energy. 

Yet another representation of the Bessel electron wave- 
function was given by Bliokh and co-workers in Ref. 0 - 
For the sake of shortness, we will not present here this— 
rather lengthy—expression. We just mention that in¬ 
stead of our Eq. JS3, constructed from the (plane-wave) 
components 0 0 with a definite helicity A, Ref. fl5l | 
employs the plane-waves whose polarization state is de¬ 
fined with respect to some (overall) quantization axis in 
the electron rest frame: 

/ y/e ± (e 2 ) \ 

Mr) = e ipr _ • (A10) 

\ V £ ~ m e (o’w) w (A )(e 2 ) / 

Although the application of Eq. (IA10I) is useful to explore 
the spin and orbital angular momentum currents in free- 
propagating vortex beams, it is less convenient for the 
analysis of high-energy electron scattering. 
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